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Abstract 

The finite size effects and relativistic corrections in pionic and kaonic hydrogen are evaluated by 
generalizing the Breit equation for a spin-0 - spin- 1/2 amplitude with the inclusion of the hadron 
electromagnetic form factors. The agreement of the relativistic corrections to the energies of the 
mesonic atoms with other methods used to evaluate them is not exact, but reasonably good. The 
precision values of the energy shifts due to the strong interaction, extracted from data, are however 
subject to the hadronic form factor uncertainties. 

PACS numbers: 13.40.Ks, 13.40.Gp 



1 



More than fifty years after their first appearance hadronic atoms continue to be im- 
portant for a better understanding of fundamental interactions. One of the first speculations 
of their existence came from the historic papers of Fermi, Teller and Wheeler 2J , where they 
showed that the time required for a negative meson (they were actually referring to muons) 
to be trapped into an atomic orbit would be (~ 10~ 13 s) much less than its mean weak 
decay lifetime (~ 10 -6 s). The negative hadron which is generally trapped into an excited 
state, undergoes transitions to lower states until it eventually enters the field of the nuclear 
strong interaction. The energy levels and widths of the hadronic atomic states are natu- 
rally affected by the strong interaction and hence experimental programmes to measure the 
shifts in the energies and widths in pionic 2I, 4], kaonic 5], S-hyperonic, antiprotonic and 
pionium [6| atoms accurately are being carried out vigorously with the aim of pinning down 
the strong interaction parameters. However, the extraction of these parameters to a good 
accuracy, requires the determination of the electromagnetic corrections accurately too. For 
example, the availability of precision data on pionic hydrogen [3[ and deuterium |7| has led 
to calculations of various electromagnetic corrections to the hadronic scattering lengths to 
better than 1% 

Whilst most of the calculations in recent literature aim at a high accuracy in evaluating 
corrections such as those due to vacuum polarization, relativistic recoil and other higher 
order corrections, the finite size of the pion and the proton is treated in a rather simplistic 
way. The correction to the binding energy of the pionic hydrogen, due to the extended 
charge of the pion and the proton is given in some works as 101 ]. 
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where, fi is the reduced mass of the up system, r n and r p the charge radii of the n and 
p respectively and a the usual fine structure constant. In sj], the Coulomb potential was 
modified by introducing a Gaussian charge distribution which depended on the pion and 
proton charge radii. In the present work, we evaluate the relativistic and finite size cor- 
rections (FSC) by modifying the Breit equation [111] to include the meson (pion or kaon) 
and the proton electromagnetic form factors. For similar Breit-like approaches, see 12]. 



The results o: 



this Breit type equation approach are compared with an 'improved Coulomb 



potential' [13] which has been used in {3, |8|, to obtain corrections due to relativistic recoil 



and the anomalous magnetic moment of the proton, while extracting the strong energy shift 
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in pionic hydrogen. Using the available parameterizations of the hadron form factors, we 
also investigate the uncertainty in the estimate of the FSC. Considering the high precision 
with which the strong energy shifts and widths for the Is pionic hydrogen states, namely, 
e u = -7.108 ± 0.013 (stat) ±0.034 (syst) eV and T ls = 0.868 ± 0.040 (stat) ±0.038 (syst) 
eV j^l, as well as the hadronic irN scattering length, a^- p = 0.0870(5) m" 1 [9] are being 
quoted and the accuracy with which the one loop calculations for the ground state energy 



of the pionic hydrogen are carried out 



3, 



the present results become relevant. 



There is no unique approach to calculate relativistic corrections to level shifts of bound 



two-body systems [12, ll3|, ll5] • We shall employ the technique of the Breit equation as it 
is particularly suited to include form-factors effects in a rather transparent way. This way 
an equation emerges which combines relativistic and finite size (FSC) effects. A further 
motivation to use the Breit approach is to compare it with results obtained in a different 
way. Regarding the relativistic corrections, it is known that the Breit equation is consistent 



at the order a 



and u sing first order time-independent perturbation theory to calcu- 



late the energy corrections jlTj] - The presence of negative energy states poses a problem in 
perturbation theory at higher orders jl7j ]. A detailed comparison between the results ob- 
tained in the Breit framework and an equation which correctly projects the positive energies 
has been performed in [Tsl. The correction to the Breit energy in this work is given as, 
AE CC = -2a 5 /j, 3 /3ttM w M p (p is the reduced mass, M n the pion and M p the nucleon mass) 
which applied to pionic atoms gives 6 x 10~ 5 eV. This is too small to be of relevance here. 

To evaluate the complete electromagnetic potential, we expand the amplitude for np 
elastic scattering, in 1/c 2 terms, thereby generalizing the Breit type equation [ll| by the 

inclusion of the proton and pion electromagnetic form factors. This leads to non-local 

□ 

terms in the potential, whose contributions are not negligible [19j. The pjp and the n ^tx 
vertices can be written in terms of the form factors F\ , F% (representing the charge and 
magnetization distributions in the proton) and F n (charge distribution in the pion) as, 

K = F'(l 2 ) (P2 + Ptf . 

The photon four-momentum, q = P[ — P\ = P2 — P' 2 - In the non-relativistic limit (q° = 0) 
and q 2 = —Q 2 , where, Q = p{ —p\ = f>2 —p{- The amplitude for the process n~ +p — > n~ +p 
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FIG. 1: Feynman diagram for pion-proton scattering 



is given by [19], 



Mji = -a[(u(ft0r^tt03i)) r£] {i/pE 2 / c ) (i/^/c), (3) 

— * 

where, D^ U (Q) is the photon propagator and u(pi), u(p{), the Dirac spinors given as, u = 



v / 2M[(l 



8Af2c 2 



), 2T?-w)] T - Substituting for the tt(p)'s and the vertex factors, and 



2M D c ' " — ~" w " — - p 

T^, the amplitude in ([3]) is evaluated and then rearranged to be written in the form, 



M fi = -2 M p [< U(pi,P2, Q)wx], 



(4) 



thus obtaining the potential in momentum space: U(pi,p2,Q) = —^it a J2i Ui(pi,p2,Q)- 
where, 



U* 
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[i ai • (Q x pi)] , £/ 9 



2 M p c 2 q 2 
Ff(Q 2 )F-(Q 2 



[i (ffi x Q) • p 2 ] 



i (pi x Bi) ■ Q] . 
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The potential in r-space is got by Fourier transforming each of the above terms 



^(Pi,p 2 ,r) = J e^U^p^Q) 



d 3 Q 
(2vr) 2 



111 ], namely, 
(6) 



The vectors pi and p 2 become differential operators in r-space ljj. The FSC to the Is state 
in pionic hydrogen can now be calculated as, AE — ^ + J2i AEi, where, 



AEi = J Wioo(r) ^(pi,p 2 ,r) ^ 100 (r)dr 



(7) 
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and ^/ioo( r ) = e~ r//ct / (a/ it a 3 ) with the Bohr radius a = 1/aji and /x the 7rp reduced mass. 
The factor a/a in AS, arises from the fact that the potential Vi(pi,p 2 , r) contains the usual 
(1/r) Coulomb potential too which must be subtracted while calculating AEi. The spin- 
dependent terms (Uq to Ug) do not contribute to AE for the Is state (see the appendix of 



19\). Expressing F p and F 2 , in terms of the Sachs form factors [l9j], G P E {Q 2 ) and G P M (Q 2 



and using Eqs (0 [6] and [7j), the total AF is given as a sum of the terms: 

fa A „\ Arn 32 a f°° F n (Q 2 ) G P F (Q 2 ) dQ a 

AE 12 = (- + AE 1 )+AE 2 = / f L )SJ + ~ 8 

\a J 7r a Jo {a' z + Q 1 ) 1 a 



tan~ 1 (Qa/2) a' 



Q a' 2 + Q 2 



dQ 



AE 3t = A£ 3 + AS = r « 

7T M p M,,- a JO 

AE ^ 7%* r WTW dQ • 

with, a' — 2/ a. In the above, the individual contribution due to AE 2 is found to be much 
smaller than (a /a) + AEi and in fact the two expressions can be combined to be written 
in the above compact form for AF 12 . The term C/3 in (jSJ) gives rise to the correction AE 3 
which cancels exactly with part of the term arising from AE4 and hence we write the total 
sum of these two terms above as AF34. Putting F p = G P E = F n = 1 in (JBJ), i.e. in the case 
of point hadrons, one gets from (jSJ) the Coulomb term plus relativistic corrections. 

We use two forms for the form factors of the proton. In the standard dipole form, 
G P E (Q 2 ) ~ G p M (Q 2 )/fi p * G D (Q 2 ), with, G D (Q 2 ) = 1/(1 + Q 2 /m 2 ) 2 , m 2 = 0.71 GeV 2 
and Hp the magnetic moment of the proton. The other parameterization is one of the 
latest phenomenological fit [^(J, where, G P E and G P M are given by the ansatz, Gm(Q 2 ) = 
G S (Q 2 ) + a b Q 2 G b (Q 2 ). The explicit forms of G S (Q 2 ) and G b (Q 2 ) are given in [20J and 
;he parameters for the proton form factors are given in Table II of 0|- The existing data 



21 



22J on the pion form factor is well reproduced by a monopole form, namely, 



1 A 2 
F n (0 2 ) = - = — (9) 

W) l + ((r 2 )/6)Q 2 Al + Q 2 {) 

such that, A 2 = 6/(r 2 ). The corrections (jHJ) can be evaluated analytically, using the dipole 
form of the proton form factors. Since the analytic expressions for AE are lengthy, we give 
below, only the leading terms (in a) of each of these terms. The sum of the corrections 
using Ui and U 2 is denoted as AE[ e 2 ad , that coming from C/3 and C/4 in (jSJ) as AE l ^ d and one 
arising due to C/ 5 as AEf ea,t 



AE- 



5 

16 am 4 A 2 
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2M 2 d 1 + 3M 2 d 2 _ A 2 n d : 



mi3 m 2 3 4m 2 3 
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TABLE I: Corrections AE in eV to pionic hydrogen using (r 2 ) = 0.439 fm 2 . Numbers in brackets 
correspond to (r 2 ) = 0.5476 fm 2 . The errors bars are due to the errors on the proton form factors. 





Ff, Ff of Ref.[20j 


Dipole form 


AE 12 (eV) 
AE U (eV) 
AE 5 (eV) 
Total AE (eV) 


0.102±0.009 (0.111 ±0.009) 
0.0388 (0.0388) 
0.0029 (0.0029) 

0.143±0.009 (0.153±0.009) 


0.095 (0.104) 
0.0388 (0.0388) 
0.0029 (0.0029) 

0.137 (0.146) 
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where we denote, m^- 



mj, with mi 



,'2 



m 2 



m 42 m 43 



4M 2 m 3 = A 2 and m 4 



(11) 
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Recall that a = 1/a/i, a' = 2/a and hence each of the above AE terms are proportional to 
a 4 . Note that if one further expands the coefficients d\ and d 2 to retain only the leading 
terms, one indeed recovers Eq. (CQ) from AE[ e 2 ad above. The calculations using the recent 
parameterization of are performed numerically. In Table I, we list the corrections to the 
binding energy of pionic hydrogen, AE, using the two parameterizations of the proton form 
factors as well as two different values of A^ in the pion form factor. The value of (r 2 ) = 0.439 
fm 2 is obtained from older we scattering experiments [21] and (r 2 ) = 0.5476 fm 2 is taken from 
a recent measurement at the Mainz Microton facility [23j . Although it is usually agreed that 
the true pion charged radius (r 2 ) = 0.439 fm 2 , we have displayed the sensitivity of the energy 
correction to the pion radius by invoking the result of a second independent measurement 
23| . As noted in 23J], the disagreement between the two measurements is supposedly due to 
a model dependence in the extraction of the value of the radius. The error on the value of 
AE 12 is evaluated using standard error propagation methods (see below). The contributions 
AE U and AE 5 are not sensitive to these errors. 



6 



In what follows, we shall compare the relativistic corrections of the present approach with 
approaches in literature which have been used for the extraction of the strong interaction 
shift, ei s , in pionic hydrogen, defined as, e is = E^; ls — E^ e Jlff ed . E^^f ed is the measured 
transition energy |3j and E^ p "^ ls is the calculated electromagnetic transition energy (here the 
strong interaction shift of the 3p state is assumed to be negligible). E^ p n ^ ls consists of the 
energy due to the Coulomb potential between point particles and various electromagnetic 
corrections [8] . Let us first consider the relativistic correction to the standard non- relativistic 
Schrodinger equation. The Hamiltonian of the Breit equation (in the centre of mass system, 
where p\= pi = — «V = p) is given as, HBreit = P 2 /2/i — p 4 /8u?c 2 + V(p,r). Evaluation 
of the second term in the above equation, treating it as usual 25j| as a perturbation, leads 
to the relativistic correction to the Bohr energy (Ejf = —pa 2 /2) of the Is state, namely, 
AE 1 ^ = — (5/8) pa 4 = — 0.215 eV. In the 'Improved Coulomb Potential' (ICP) approach 
of Ref. [I3j, starting from Eqs (9) and (10) in 13|. one can find the total energy, Ejf + E ICP , 
for the case of a spin-1/2 and spin-0 bound state, where, 



E 



ICP 



pa + 2n r 



p y f pa 



Mj V 2 



pa 



p 



2p 



- 0.215 eV + 0.01 eV + 0.037 eV 



4(M p + M n ) M p + M n 

.Is 

rel 



2 n 



AE 1 * + 0.047 eV. 



(12) 



Here, k p = p p — 1, with p v = 2.793 nm. The first term represents the relativistic correction, 
which is referred to as the standard Klein-Gordon result in jl^ . This term is identical to 
AElli obtained from the Breit type equation. In the third approach, one could actually use 
the Klein-Gordon (KG) equation as was done in 3j, |8j. Here the difference between the KG 
result, E l K s G and the Bohr energy, Ejf is, E l K s G - Ejf = AE l K s G = -0.211 eV. 

In order to compare the terms apart from AE re \ in the Breit type equation approach with 
those in [l3j, we assume point-like hadrons such that the energies in (jSJ) become, 



AE\l + AE ls 



a 4 M p Ml 



1 + 2^ 



0.0417eV. 



(13) 



2 (M p + M^f 

with the tilde indicating the fact that the energies correspond to point-like hadrons. From 
Table I, we can see that AE 3i + AE 5 is not different from A£ , 34 + AE 5 (up to the fourth 
digit after the decimal) and the effect of the hadron form factors on these two corrections 
is negligible. Besides this, we also note that in contrast to 131 ] . the contribution of the 
proton magnetic moment in the present work is found to be negligible. This can be seen by 
examining Eqs fflOl) which are obtained analytically assuming dipole proton form factors. 



TABLE II: Contributions in (eV) to E^i^ ls , and the deduced strong interaction shift, t\ s usin^ 
(r 2 ) = 0.439 fm 2 . Numbers in brackets correspond to ((r 2 ) = 0.5476 fm 2 ). 



Point Coulomb, E 3 B p ~* ls + AE%~ 


► Is 


2875.7196 


Breit type equation (with finite size) 


-0.143±0.009 (-0.153±0.009) 


Vacuum Polarization, order a 2 [I 




3.235±0.001 


Vacuum Polarization, order a 3 


[3J 


0.018 


Vertex correction [3j 




-0.007±0.003 


Pionic atom recoil energy [3j 




-0.004 


Total calculated E^^ ls 




2878.8186 ±0.009 (2878.8086 ±0.009) 


^measured yjj 




2885.916 ±0.013(stat) ± 0.033(syst) 


, jpe.m. rpmeasured 




-7.097 ±0.009(FSC) ±0.013(stat) ± 0.033(syst) 
(-7.107 ±0.009(FSC) ± 0.013(stat) ± 0.033(syst)) 



Though, Ff in (jSJ) contains both the electric and magnetic Sachs form factors, there appears 
no term with k p in the corrections at leading order in a as in ffTU]) . To summarize the above, 
we have three different approaches of summing the relativistic corrections: 



AE^r ls - 0.0417 eV 



E%~; ls = AE 3 r y ls - 0.047 



j-i3p^ls 
^Sigg 



AE^ U 



0.047 



0.171 eV 
0.166 eV 
0.161 eV. 



(14) 



As can be seen there is a slight dependence on the approach used to calculate the relativistic 
corrections. It is somewhat inconsistent to use E Sl " j^, 8] as the sum of relativistic correc- 
tions, since E Sl " is a sum of AE^^ ls and 0.047 eV, where 0.047 eV is taken from E ICP 
(where AE^![* * ^ AE^q 1s ). Using a correction of AE^ t = 0.0239 to the Bohr energy 
of the 3p state, namely, E^? = 359.441, in Table II we present a consistent deduction of 
the strong energy shift. The relativistic and FSC are taken from the Breit type equation 
approach and the remaining corrections are as in With the potential ([6]) being short- 
ranged, the finite size and relativistic corrections, AEi 2 + Ai? 34 + AE 5 , to the energy of the 
3p state are very small and hence neglected. 

As evident from Table II, the error due to the electromagnetic form-factors of the proton 
is of the same order as the statistical and systematic counterparts. Therefore some remarks 
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on its determination are in order. The 6x6 correlation matrix pij = Cov(a,i,aj)/(Ji(Jj (a.j 
are the fitted parameters and cr, their respective errors) was supplied to us by the authors 



of 



20l | . The error on AE due to uncertainties of hadronic form-factors is calculated by the 



standard method, i.e. 

(15) 



M 1 



uovfaj, a,-)— — 
o oaj 



where the subscript denotes the central value. Taking (A%) 2 = 1, we obtain the l-er error 
on AE, namely, (5E)p a sc = ±0.009 eV. For 2-cr variations in the parameters, A% 2 increases 
by 4 and the error on AE is doubled. 

Within the framework of the present work, the correction to the energy of the Is state 
in kaonic hydrogen (using the proton form factors of [20] and a monopole kaon form factor 
with (r 2 K ) = 0.34 fm 2 ) is, E^ Q f = + AE (FSC) = 0.573 eV + 2.525 eV = 3.098 eV 

(using central values of form factor parameters). This correction would be relevant when 
better data on kaonic hydrogen would become available from the ongoing programme of the 



DEAR collaboration [5]. 

In summary, we can say that the present work investigates the relativistic and finite size 
corrections in hadronic atoms, using a Breit-type equation. These corrections have been 
shown in the present work to be important for the precision measurements of the strong 
energy shifts in pionic hydrogen. We find that the contribution of the magnetic moment of 
the proton to the corrections is negligible. In future, we plan to extend such calculations 
for the evaluation of a spin-0 - spin-1 amplitude which would be relevant for the pionic 
deuterium case. The full electromagnetic potential in the ird case will involve the deuteron 
electric, magnetic and quadrupole form factors. In the ird atom, the strong energy shift has 
been found to be repulsive, namely, ei s = 2.43 ± 0.1 eV [7], with the contribution of the 
FSC, 0.51 eV (using the simple Eq. ([1]) with the proton radius replaced by the deuteron 
radius). The above approach could alter the precision values for pionic deuterium obtained 
so far. 
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